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5 We discuss several constraints on the properties of 7-ray bursts (GRB) at cosmological 
distances. First we use the requirement that burst sources must be optically thin to a 
test photon for the process 7 + 7 — > e + + e~ in order to produce the observed nonthermal 
spectra. In particular, we derive probability distributions for the minimum Lorentz expansion 
factor 7 m i n , the radiation energy E~ /} the maximum baryonic mass Af max , and the maximum 
surrounding gas density n max in the events, based on 254 events from the second BATSE 
catalog. In the case where the GRB spectrum cuts off at the highest observed energies 
(~ 100 MeV), we obtain the mean values (7 m i n ) = 90, (E 7 ) = 4 X 10 51 7i — 2 ergs, and 
(M max ) = 3 X 10 -5 <f; -1 M©, where £ is the fraction of the total energy which is converted 
to 7-rays. The distribution of burst energies ends at about 10 53 ergs, close to the binding 
energy of a neutron star. 

Secondly, the time variabilities of the bursts in the BATSE catalog are used to place an 
upper bound i? max on the curvature radius of the emitting surfaces in the events. This is 
based on the requirement that the emitting region seen by the observer must be sufficiently 
small to produce the observed variability without violating causality. Using the 64-msec 
resolution of BATSE, we find that a significant number of bursts have i? max ~ 10 13 (7/10 2 ) 2 
cm, where 7 is the Lorentz factor of the expansion. This limit should become stricter with 
finer time resolution. 

Finally, we discuss the association of cosmological GRB with galaxies. We consider 
eight bright and well-localized bursts detected by the Pioneer Venus Orbiter, whose posi- 
tional error boxes contain no bright galaxies. Assuming that burst events occur in galaxies, 
we place upper limits on the luminosities of the host galaxies. Using the local luminosity 
function of galaxies, we calculate the probability for not seeing the GRB host galaxy. This 
probability tends to increase as the width of the GRB luminosity function increases. How- 
ever, the allowed width of the GRB luminosity function is restricted by the burst peak flux 
distribution. 

Subject headings: cosmology: observations-gamma rays: bursts 
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1. INTRODUCTION 

Gamma-ray bursts (GRB) have eluded a definitive explanation for more than twenty 
years. Hundreds of models have been proposed in the past, most of these involving Galac- 
tic sources (cf. Nemiroff 1994). More recently, however, the Burst And Transient Source 
Experiment (BATSE) on board the Compton Gamma-Ray Observatory (GRO) has shown 
that the burst population is highly isotropic (Meegan et al. 1993; Briggs et al. 1993). This 
result is inconsistent with a Galactic disk population of sources, but suggests that the bursts 
occur either in an extended Galactic halo or at cosmological distances (Mao & Paczyhski 
1992; Hakkila et al. 1994). There is also a deficiency of faint bursts relative to a uniform 
distribution in Euclidean space, which can be naturally explained if the faint bursts are at 
high redshifts (Mao & Paczyhski 1992; Piran 1992; Dermer 1992). A cosmological redshift 
is also implied by the suggestion that dim bursts are on average longer than bright bursts 
(Norris et al. 1993). Since little is known about the origin of cosmological GRB at this stage, 
it is important to extract the maximum amount of information on the burst sources from the 
existing data. In this paper we discuss general empirical constraints on cosmological 7-ray 
burst models, and use the second BATSE catalog to quantify these constraints for the entire 
burst population. Some of our constraints were briefly discussed in a previous paper on the 
subject (Woods k Loeb 1994). 

First, let us consider qualitatively the physics behind the various constraints derived in 
this work. A typical burst located at a distance of order the Hubble radius (c/Hq ~ 3/i -1 
Gpc) must release a large amount of energy (~ 10 51 ergs) to produce the observed fluence. 
Combining this with the small source size (^ 3 X 10 10 cm) associated with the typical variabil- 
ity timescaleof ~ 1 second, one obtains an optical depth to pair production (7+7 — > e + -\-e~) 
of order ~ 10 11 . Under these circumstances, the photon distribution should be thermalized, 
in conflict with the observed nonthermal spectra of GRB (vl v const between 0.1 and 100 
MeV). The thermalization problem is avoided if the the emitting photosphere of the burst 
source expands relativistically with a Lorentz factor 7 ^> 1 (Goodman 1986; Paczyhski 1986). 
The emission is then beamed into cones of half angle ~ 7 -1 , so that the center-of-mass mo- 
mentum in photon-photon collisions is lowered below the pair-production threshold, and the 
characteristic photon paths cross only at large distances where the scattering probability is 
low. For typical burst parameters, the requirement for an optically thin source translates to 
the bound 7 <; 10 2 for emission from a fixed radius (Fenimore, Epstein, & Ho 1993). For 
emission from an expanding shell, this constraint is weakened considerably, so that 7 ^ 10 2 
is allowed. 

This limit on the Lorentz factor may be used to place constraints on other burst prop- 
erties. The minimum 7 can be used to set an upper bound on the baryonic mass M partic- 
ipating in the GRB event (Shemi & Piran 1990). This constraint is obtained through the 
relation M = Ej/^jc 2 , where Ej is the total 7-ray energy released, c is the speed of light, 
and £ is the fraction of the total energy released in 7-rays. Based on this upper bound on 
the number of baryons in the burst, it is possible to limit the ambient gas density which is 
swept by the expanding photosphere of the source. We define this bound as n max . 

Until recently all of the above limits were evaluated only roughly for "average" values of 
the burst parameters. In an earlier publication, the first BATSE catalog allowed us to get 
a first statistical evaluation of these constraints for the entire burst population, assuming 
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emission from a fixed radius (Woods & Loeb 1994). In §2.1 and §3.1 we assume the more 
physical case of emission from an expanding shell, and derive these constraints for an even 
larger sample of bursts, using the second BATSE catalog and its associated data set of GRB 
light curves. Moreover, we include two additional constraints as described below. 

Some models require that the main part of the emission take place at a radius larger 
than some minimum radius. In the fireball model of Meszaros, Laguna, & Rees (1993), a 
burst is produced when the expanding fireball decelerates due to the surrounding medium. 
For typical burst energies (E ~ 10 51 erg) and typical interstellar densities (~ 1 cm -3 ), the 
fireball is not decelerated appreciably until it has reached a radius R ~ lO^^/lO 3 ) -2 / 3 cm. 
If an empirical upper bound can be placed on the radius at which the burst is produced, 
then such models can be constrained. 

An upper bound of this type can be obtained from the variability timescales of the 
bursts in the second BATSE catalog. Let us assume that the emission takes place at a 
radius R while the emitting surface expands with speed /3c and corresponding Lorentz factor 
7 = (1 — /3 2 ) -1 / 2 , and that signals propagate in the rest frame of the expanding material at 
speeds ~ c sg . A burst which exhibits variability SC/C in its count rate C on a comoving 
timescale St' must have originated from a region smaller than the signal-travel time £ ~ c S gSt' 
due to causality. The observed timescale corresponding to St' is St = j(l — /3)St' (2j)~ 1 St l , 
where the first factor 7 is due to the Lorentz transformation from the comoving frame to the 
lab frame, and (1 — (3) accounts for the velocity difference between the expanding material and 
the photons. Hence, £ ~ 2jc S gSt . An area larger than ~ £ 2 will in general contain many 
regions which are causally disconnected from one another on the observed timescale St . 
These regions act as independent emission sites, and by Poisson statistics produce variations 
SC/C ~ A^ -1 / 2 , where N is the number of emitting regions of size ^ £ that are effectively 
contributing to the observed signal. Since relativistic beaming limits the emission to an angle 
~ I/7, the effective emitting area of the source is ~ 2 R 2 , where R is the curvature radius of 
the (not necessarily spherical) emitting surface. Therefore, the minimum number of causally 
disconnected regions at the emission photosphere is given by N ^ 6 2 R 2 / £ 2 ~ (R/2j 2 c S gSt ) 2 , 
and the resulting fluctuation amplitude is SC/C ^ 2j 2 c S gSt / R. This condition can in 
turn be used to find the maximum curvature radius i? max = 2j 2 c S gSt / (SC / C) which could 
produce the observed variabilities SC/C in each individual burst. We discuss this constraint 
in detail in §2.2 and §3.2. 

As the last empirical constraint, we consider the potential association of GRB with 
galaxies. Searches for GRB counterparts at other wavelengths have all turned up empty- 
handed (Schaefer 1993), and there is doubt as to whether extragalactic bursts could even 
be associated with galaxies (Schaefer 1992; Fenimore et al. 1993). Any model which places 
burst sources inside of galaxies must be consistent with observational constraints on the 
galaxies near GRB positions on the sky. If bursts occur outside of galaxies, one needs to 
turn to exotic emission mechanisms (e.g. cosmic strings; Paczyhski 1988; Babul, Paczyhski, 
& Spergel 1987) or postulate the existence of an as yet unobserved population of neutron 
stars in the intergalactic space. In the latter case, the neutron stars could have been ejected 
from galaxies with large kick velocities (Narayan et al. 1992). 

The positional error boxes associated with eight well-localized bursts were searched for 
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galaxies, and none were found down to very faint magnitudes B ^ 24. This null result was 
used by Schaefer (1992) to derive a lower bound on the distance to a burst, assuming the 
burst is associated with a galaxy whose luminosity is comparable to that of M31. Fenimore 
et al. (1993) used Schaefer's results, together with the combined BATSE-PVO peak flux 
distribution, to obtain limits on the absolute magnitudes of the GRB host galaxies. In 
§4, we derive more restrictive constraints from the same data. In particular, the galaxy 
luminosity function is used to set constraints on the width of the GRB luminosity function. 
Obviously, if the luminosity function of GRB is sufficiently broad, it should allow for some 
bright bursts to be at large distances, where their host galaxies are too faint to be seen. 
Thus, the assumption that most GRB occur inside galaxies can be used to set a lower limit 
on the width of the GRB luminosity function. 

A summary and discussion of all of the above results is given in §5. Throughout the 
paper we assume a Hubble constant Hq = 75 km sec -1 Mpc -1 . 

2. CONSTRAINTS ON SOURCE PROPERTIES: THEORY 

In this section we derive constraints on properties of the burst sources in terms of quan- 
tities measured by the BATSE experiment. The distances to the individual bursts are cali- 
brated according to their peak flux. This distance estimator provides an excellent fit to the 
number count statistics in standard cosmologies (cf. Fig. 4a, later). 

2.1 Pair Production Optical Depth of Cosmological Bursts 

We first calculate the optical depth for pair production by photon-photon collisions near 
a cosmological source. Because of the conservation of relativistic energy-momentum, a test 
photon of energy et can only produce an e + e~ pair in a collision with a photon whose energy 
e is greater than eth = 2mgC 4 /et(l — cos#), where m e is the electron rest mass and is the 
angle between the photon trajectories. The cross section for collisions above this threshold 
is 



( 3 _ v ± ) lJA±^l_2v(2-v 2 ) 



where <r T is the Thomson cross section and v = [1 — (e t h/e)] 1//2 is the center-of-mass speed 
of the outgoing pair particles in units of the speed of light (e.g. Berestetskii, Lifshitz, & 
Pitaevskii 1982). 

Suppose now that a large amount of energy is released at a point in space, in the form 
of relativistically expanding matter (baryons, e + e~ pairs, or even exotic particles). At some 
distance from the source, the matter begins to produce 7- rays isotropically in its rest frame; 
the emission continues until the expanding shell has reached a radius i?o (see Fig. 1). The 
mechanism of 7-ray production is not important; it could be the interaction of a baryon- 
pair fireball with the surrounding medium (Meszaros et al. 1993 and references therein), 
or collisions between shells (Rees & Meszaros 1994). If the matter expands with a Lorentz 
factor 7^1, the emission is strongly beamed, and the observer will see photons only from a 
small patch of the emitting surface. Thus, for our purposes the emitting surface need not be 
spherical; it could be a collimated jet with an opening angle > 7 -1 and a curvature radius 
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Fig. 1: Geometry for the calculation of the pair-production optical depth 
r for an expanding source. The shaded area represents the contribution to 
r at a distance R from the source from a ring of angle ip and angular width 
dip. The surface shown is the locus of points from which photons originate 
whose trajectories intersect that of the test photon at the instant it reaches 
a distance R (Fenimore et al. 1993). 



Consider a test photon of energy et which is emitted at the end of the emission, at a 
distance i?o from the center of the explosion, and which then propagates precisely along the 
line of sight. When the test photon has reached a distance R from the center, only photons 
which arrive at this point at precisely the same instant can collide with the test photon to 
produce pairs. The locus of points from which these photons arrive defines a surface i? e (?/>), 
where R e is the emission radius and ip is the angle the emission radius vector makes with 
the line of sight at the center of the explosion (see Fig. 1). This surface is given by 

(R 2 e -2RR e cos^ + i? 2 ) 1 / 2 , -Re _ -Ro , R-Ro (e) , 

+ /3c ~ /3c c ' { ) 

where (3 = (1 — 7 -2 ) 1 / 2 to be the expansion speed in units of the speed of light. The solution 
to equation (2) is 

R e = 7 2 [i? + P(R-Rq) - /3 2 i?cosV>] 

- { 7 4 [i? + P(R - Rq) - /3 2 i?cos V>] 2 - 7 2 (1 " PW - P)Rq + 2/3i?i? ]} 



1/2 (3) 



We assume that the emission is steady over a timescale ~ Ro/c. The optical depth r 
is obtained by integrating the collision probability from i?o out to infinity. At each point 
in the test photon's path, one must also integrate over the photon energy distribution, and 
over the contributions from different collision angles: 



oo oo 



r(e t ) = ^JdRjdej dA **^ a{e t ,e,0), (4) 
where d<p(e, R) is the contribution to the photon flux per unit energy, at a radius i?, from 
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an area element dA on the emitting surface. Photons emitted from dA collide with the test 
photon at an angle 6. 

Figure 1 shows the contribution at a distance R to the optical depth from photons 
originating in a ring which subtends angles in the range (ip,ip + dip) at the center of the 
explosion. We define 8 to be the angle between the photon trajectory and the radial direction 
in the observer's frame. Note that the observed energies of photons from this ring are blue- 
shifted from their rest-frame energies by a factor [7(1 — /3cos <5)] _1 . However, for 7 ^> 1, the 
surface "seen" by the test photon at any given instant always covers an angle ~ 7 -1 , due to 
the restriction imposed by equation (2), and the effect of relativistic beaming (cf. Fenimore 
et al. 1993). Thus, the range of blueshifts contributing to the total flux is the same at all 
distances R from the source, so we may make the approximation that the same spectrum 
is observed at all radii. A detailed calculation shows that this approximation yields nearly 
the same results as the case where the variation of the spectrum with radius is taken into 
account. 

The photon flux per unit emitting area at a distance R from the source is given by 

d<f>(e,R) <f>(e,R) 



dA J(R) [7(1 - /3 cos 8)\< 



(5) 



The two factors of 7(1 — /3cos<5) account for the relativistic beaming due to the Lorentz 
transformation of solid angles (cf. Rybicki & Lightman 1979). We have defined the normal- 
ization J(R) such that the integral of d(p/dA over the entire contributing surface yields the 
total flux </>(e, R), 

T(R) _ f dA f 2irR 2 e [l + (R' e /R e ) 2 ] 1 / 2 sm^d^ 

1 } ~ J [7(1-/3 cos 8)] 2 J [7(1-/3 cos 8)] 2 ' 1 ] 



where R e is given by equation (3), R' e = dR e /dtp } and ip ma , x (R) is the maximum angle which 
contributes at a distance i?, i.e., the angle such that the emitted photon's trajectory is 
tangent to the surface (see Fig. 1). The flux per unit energy at a distance R is given in 
terms of the observed flux per unit energy <^ (e) by 

^ R) = ^r 7T^—^. to (rrz) > (7) 



R 2 (1 + z) 2 Y ° \l + z, 

for a source at a cosmological redshift z, and a corresponding luminosity distance D(z) (e.g. 
Weinberg 1972). The additional factor of (1 + z)~ 2 is due to the fact that we are considering 
the number flux of photons per unit energy rather than their total energy flux. Equations 
(l)-(7) then provide the total optical depth of a GRB source, with a bulk expansion Lorentz 
factor 7, to a test photon with energy et: 

00 00 

T (7> e t) = — — — 77 / D o J( m / de 



c{l + z) 2 J R 2 J{R) J ru \l + z 



^max . 

. 2TrR 2 e [l + (R l jR e ) 2 ] 1 / 2 
# Sm ^ [7(1-/3 cos £)P ^'^ 
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Note that the cross-section a depends on the collision angle 6, while the relativistic emission 
probability depends on the emission angle 6. Inspection of the triangle in figure 1 yields the 
relations 5 = ip -\- 6, tan 6 = Ro sin ip/(R — Ro cos ?/>), and sin 6 = (Ro / R) sin 6. 

The burst spectrum can usually be approximated by a power law, 

<f) (e) = kc~ a } ei < e < e 2 , (9) 

Most burst spectra seem to be more or less featureless continua with roughly equal energy 
flux per logarithmic energy interval, i.e., a = 2 in equation (9). The spectra are highly 
nonthermal, with no sign of a significant decline at any value of £2, even out as far as 
£2 = 100 MeV. For an expanding source with a = 2, the optical depth r is quite sensitive 
to the energy cutoff e2. Fenimore et al. (1993) showed that for large enough values of 7, 
all of the photons will collide with the test photon at angles less than the critical angle 
9 CV 2m e c 2 /(ete2) 1//2 , meaning that the center-of-mass energy will always be below the 
threshold value eth and no pairs will be produced. Since burst spectra have been observed to 
extend out to energies of order 100 MeV, we take a test photon energy et = 100(1 + z) MeV, 
and will consider two extreme cases: the case where the spectrum extends out to infinite 
energies (e2 — > 00), and the case where the spectrum cuts off at an energy of the order of 
the test photon energy (e2 = 100 MeV). Note that since most of the contribution to r comes 
from near the source, we are justified in integrating out to infinite radius. 

Implicitly built into this calculation is the assumption that the emission is steady on 
timescales ~ Ro/c. This corresponds to an observed timescale 8t ~ (1 — /3)Rq/c 
i?o/(2c7 2 ), where the factor (1 — f3) accounts for the shortening of the duration due to 
the source's expansion toward the observer. This is of the same order as the geometric time 
delay between the arrival of photons along the line of sight and photons coming in from the 
edge of the observed emission region, at angles tp ~ 7 -1 . Since all signals are smeared out 
by this delay, there is no way to test this assumption empirically. However, the assumption 
is natural since plausible physical processes that convert the kinetic energy of the fireball 
to 7-rays after a time Ro/c (e.g. deceleration by the ambient gas or shell collisions) should 
typically last a time <; Ro/c. 

To satisfy the observational constraint that the burst spectra are nonthermal, we require 
r < 1. By combining equations (8) and (9) with a = 2, and imposing the inequality r < 1, 
we may obtain a lower bound Ro > -R m in(7? e t): 

where ^(7) is a numerically-calculated function of 7. In the case where the spectrum extends 
to infinite energies ( 62 — > 00), this function is very well fit by a power law, ^(7) = 0.0337 -2 ; in 
the case where e 2 = 100 MeV, it is fit to ~15% accuracy by = 0.0327- L79 [l-(7/97)] 3 - 07 . 
Figure 2 shows the numerically-calculated values of i? m i n for a burst at a distance D = 1.2 
Gpc, for the two cases. Note that for any value of 7, the 7-ray source must be in thermal 
equilibrium at radii smaller than -R m in(7)- 

One needs to relate Ro to observable quantities in order to obtain a useful constraint 
on 7 from equation (10). Such a relation is obtained from the comment made above that 
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any variation in the observed burst flux must be longer than the geometric time delay 
between photon paths (cf. §2.2). Taking into account the cosmological time dilation, T <; 
i?o(l + - 2; )/(2c7 2 ). This in turn gives an upper bound on i?o, 

2c 7 2 T 

-Rmax = -j— , (11) 

where T will be taken to be the observed duration of peak flux emission for the burst. 
Equations (10) and (11) lead to a lower bound on the Lorentz factor of the expanding 
material, 7 > 7 mm (et) = [(1 + z)Rmin/ '(2cT)] 1//2 (see Fig. 2). For the case where £2 — > 00, 
one can solve analytically for 7 m i n to obtain 



7n 



3 <r T ke t 0.033(1 + z)D 2ll/A 



32 c 2 m 2 c 4 T 



12) 



In the case where €2 = 100 MeV, equations (10) and (11) must be solved numerically for 
7min by Newton- Raphson iteration. As we will see in §3.1, the quantity kD 2 is independent 
of the observed properties of a given burst. 

Some GRB models have an unavoidable amount of baryonic matter present, such as those 
involving neutron stars (Woosley 1993; Eichler et al. 1989; Meszaros & Rees 1992; Narayan, 
Paczyhski, & Piran 1992; Usov 1992; Loeb 1993). Thus, an interesting empirical constraint 
would be an upper bound on the baryonic contamination allowed for a burst source. The 
total bulk kinetic energy of the matter is jMc 2 , where M is the total rest mass. It may be 
that not all of this energy is converted to 7-rays, and only a fraction £ = Ej/(jMc 2 ) of the 
total kinetic energy is radiated as 7-ray energy Ej. The maximum baryonic mass allowed is 
then Mmax = E 1 / (^ m [ n c 2 ). 

Also of interest is the type of environment the GRB sources inhabit. If bursts occur in 
galaxies, the density of the surrounding medium should be typical of the interstellar medium 
(~ 1 cm -3 ). The above bounds on 7 and M can be used to set an upper limit on the 
ambient gas density. For a spherical expansion, the minimum amount of matter present is 
the amount swept up by the expanding fireball, M ISM = (47r/3)mn ISM i?Q, where m is the 
mean atomic mass and n ISM is the interstellar particle number density. The requirement 
^ism < Mmax yields an upper bound on the gas density surrounding the source: n ISM < 
"max = M/(47ri? min /3), out to (Ro) £ 10 16 cm. Note that n max depends on the Lorentz 
factor 7 of the source; if GRB spectra extend to infinite energies, then to leading order 
"max oc 7 s . In §3.1, n max is calculated for 7 = 7 m i n ; however, the strong 7-dependence 
considerably weakens the constraint. 



2.2 Derivation of the Maximum Curvature Radius 
of the Emitting Surface 

Consider GRB emission from a surface which moves outward with a bulk Lorentz factor 
7 ^> 1. The emission is beamed into cones of half-angle 7 -1 <C 1. The surface need not 
be spherical; the emission could be coming from a collimated jet of opening angle <; 7 -1 or 
a part of a more complicated surface. Suppose that a particular burst exhibits significant 
intensity variations (SC/C ~ 1) in its count rate C on a timescale St' as measured in the 
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1 10 100 1000 

Lorentz factor y 



Fig. 2: Method of obtaining 7 m i n for a burst. Equation (10) is used to obtain 
i? m in(7) for (a) spectra extending to infinite energies (upper solid curve), and (b) 
truncated above 100 MeV (lower solid curve). Equation (11) is used to obtain 
-Rmax(7) (dashed curve). The intersection of i? m i n and i? max yields 7 m i n for the 
two cases. The curves shown are for a burst with an observed duration T = 1 
sec, for a test photon energy et = 100 MeV. The value of i? m i n is independent 
of the burst's peak flux (§3.1). 



rest (comoving) frame of the expanding fluid. We then postulate that the linear size of the 
emitting region for this GRB cannot be larger than c S gSt' = j~ 1 c S gSt, where St = jSt 1 is 
the time interval measured in the lab (nonexpanding) frame, and c sg is the speed at which 
signals (e.g. sound waves) propagate in the comoving frame. A region of this size is causally 
disconnected from other such regions on a lab-frame timescale St. We assume that there is 
no modulation of the source which affects all of these regions simultaneously. 

Suppose that the emitting surface has a radius of curvature R with respect to the source 
(see Fig. 3). A region of solid angle o?fi <C 1 with respect to the center of curvature then 
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Fig. 3: Geometry for the calculation of the maximum curvature radius 
i? max . After a time St, the radiation which was emitted at the beginning of 
this period has moved ahead of the matter by a distance (1 — f3)cSt. The 
causally connected regions are shaded. 



contains dN effectively independent emitting sites, where dN is given by the inequality 

dN > (^-) ■ i] dfi, (13) 
\c sg StJ 

and T) is the fraction of the surface covered by the emitting "cells" . 

The emitting region seen by a distant observer can be divided into circular rings centered 
on the line of sight. A ring of angular size ip (with respect to the center of curvature) and 
width dip contains = (dN/dft) ■ 2tt sin ip dtp emitting cells, where dN/dft is given by 
equation (13). Due to relativistic beaming, the contribution to the total mean count rate C 
from this ring is = (dC / dft) ■ 2tt sinip dip, where 

'!£ = .J £ (Mi 

dtt 47T 7 2 (1 - (3 cos iP) 2 y ' 

(cf. Rybicki & Lightman 1979); 7 ^> 1 is the bulk Lorentz factor, and (3 = (1 — 7 -2 ) 1 / 2 . For 
simplicity, the cells may simply be treated as being in either an "on" state or an "off" state 
during this time. The emission from the ring is then governed by the Poisson statistics of 
the independent causally disconnected cells, so that the variance of the intensity is (SC)^ = 
C^/N^p. The variances of all the rings add in quadrature, so the total variance is 

«* c » 2 =/(lT) 2 (f)"- 2 — < 15 > 



The quantity dN/dft can be treated as a constant (cf. Eq. (13)), and may be pulled outside 
of the integral. Combining equations (13), (14), and (15) and performing the integration 
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^ = ( 1 +-)(^2+^2-)- (17) 



yields an upper bound for the curvature radius R < i? ma x, where 

_ Cap St C , , 

R - = ^ sc- (16 » 

The lab-frame emission timescale St is related to the observed variability timescale St 

by 

St R SC 

The first term is due to the velocity difference c(l — /3) c/(2j 2 ) between the expanding 
material and the photons (cf. Fig. 3), and the second term is the geometric time delay given 
by equation (11). The factor of SC/C is necessary because only the fraction SC/C of the 
photons which contribute to the observed fluctuation suffer this geometric delay. The factor 
(1 + z) accounts for time dilation in the case of cosmological sources. In terms of observed 
quantities, equation (16) becomes 

^ max - V3^+(c sg /c) 1 + z SC' ( * } 

Notice that, to within a factor of order unity, equation (18) resembles equation (11) 
under the substitutions T ~ St / (SC /C) and c sg ~ c. However, while equation (11) was 
based on a geometric time delay, equation (18) is based on causality and could in principle 
yield stricter bounds on i? max if c sg <C c. The method for obtaining St / (SC /C) from the 
BATSE light curve data, as well as the problems involved in constraining 7, are discussed 
in the next section. 



3. CONSTRAINTS ON BURST PROPERTIES: RESULTS 

In this section we apply the constraints derived in §2 to the BATSE burst population. 
The methods of extracting the relevant information from the data are discussed, as well as 
the assumptions made about the population of burst sources. 

3.1 Distributions in E y , 7 m i n , M meLX , and n max 

The goal of §2.1 was to derive the values of 7min? ^max, and n max , and to express them 
in terms of observable quantities. We have already assumed that the burst spectra can be 
approximated by a power law with spectral index a = 2 (Eq. (9)). With this assumption 
in mind, the three quantities to be obtained from the BATSE data are the distance to the 
source D } the duration of peak emission T, and the spectral constant k. 

To obtain the distance to a burst, one must make assumptions about both the luminosity 
function of the bursts, and the cosmological model which describes the geometry of the 
universe. The GRB number count data have been shown to be consistent with a nonevolving 
population of standard candles in peak luminosity, all with identical power-law spectra <f> oc 
e -2 , in an = 1, A = universe (Mao & Paczyhski 1992; Piran 1992). In principle, one 
might also consider bursts as standard candles in total energy and use fluence as a distance 
indicator; however, because BATSE is sensitive to peak flux, this introduces a significant 



11 



selection bias against faint, long bursts. For A = 0, the statistical significance of the number- 
versus-peak-flux fit depends only very weakly on the value of (Wickramasinghe et al. 1993). 
We therefore adopt the simplest model and assume a nonevolving population of standard 
candles with a spectral index a = 2 in an = 1 universe. For a given value of a, the 
observed peak number-flux of photons C^^z) from a source at a redshift z is given by 

} r • (i + z) 2 ~ a 

Cpeak = / H e o)de = , (19) 

where e L < e < e v is the energy range picked up by the detector, Y is the intrinsic peak 
emission rate of photons per unit time of the source in this energy range, and D is its 
luminosity distance (Weinberg 1972), 

D = ^[(1 + z) - vT + 7]. (20) 

There are two kinds of peak flux data listed in the BATSE catalog: the peak flux C pea k 
in photons cm -2 sec -1 , and the ratio of maximum to threshold count rates C max /C m i n . We 
used C pea k as our distance indicator, since this quantity was corrected for detector orientation 
and atmospheric scattering. Values of C pea k are given for three integration times: 64, 256, 
and 1024 msec. The value in the 64-msec channel is likely to be the best measure of 

the true peak flux. Due to detector inefficiency at low fluxes (^ 1 cm -2 sec -1 ; Fishman et 
al. 1994), only bursts with C peak = C^ k > 1 cm 2 sec 1 were included in the sample. To 
determine the value of T in equation (19), a one-distribution Kolmogorov-Smirnov (K-S) test 
was performed on the number- versus-peak-flux distribution for the 254 bursts with C pea k > 1 
cm -2 sec -1 and no data gaps (see below). The best-fit value was 

sec -1 , where the upper and lower bounds denote the 95% confidence region, and h = i/o/(100 
km sec -1 Mpc -1 ) is the dimensionless Hubble constant. Figure 4a shows the fit obtained 
for the number counts vs. peak flux with this best-fit value of Y. Taking this value, the 
luminosity distance to a standard-candle source with peak flux C pea k is given by equation 
(19). Figure 46 shows the K-S significance for standard candles as a function of T, indicating 
that there is actually a range of values which give a good fit. It should also be noted that 
the significance of the K-S fit was insensitive to the width of the burst luminosity function, 
as already pointed out by Lubin & Wijers (1993), Ulmer & Wijers (1994), and Cohen & 
Piran (1995). A luminosity function with a width of the order of its mean is not ruled out, 
and in this case a value of C pea k would correspond not to a single value of D but rather to 
a probability distribution in values of D. Figure 4c shows the K-S significance as a function 
of "width" a for a log-normal luminosity function peaked at T = 1.0 X 10 57 /i -2 sec -1 (see 
Eq. (25)). 

The spectral constant k is obtained by assuming that the spectral shape (9), with a = 2, 
applies during the period of peak emission, and using the first equality in equation (19). The 
BATSE peak flux data are for the channels 50-300 keV, so setting e L = 50 keV and e v = 300 
keV, one has 

k = (0.06 C peak ) MeV cm -2 sec -1 . (21) 
Note that since k oc C pea k and D 2 oc (C pea k) -1 , the optical depth, and hence 7 m i n , is 
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Fig. 4: (a) Number count N vs. peak flux C pea k for the 254 bursts in the 
second BATSE catalog with C pea k > 1 cm _2 s _1 and no data gaps (solid 
line), and the best-fit curve for standard candles in an = 1, h = 0.75 
universe (dotted). (6) Significance of the K-S fit as a function of peak 
emission rate T near the best-fit value of T = 1.0 X 10 57 /i -2 sec -1 , with 
= 1 and h = 0.75. (c) Significance of the K-S fit as a function of width 
a for a log-normal luminosity function peaked at T = 1.0 X 10 57 /i — 2 sec -1 
(Eq. (25)), with = 1 and h = 0.75. 
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independent of C pea k to within a factor of order unity (cf. Eqs. (10)-(12)). This reflects 
the assumption that GRB are standard candles in peak flux, and hence all have the same 
photon density at the source. 

The peak emission duration T was formerly expressed only in terms of an upper bound, 
given by the total burst fluence (time-integrated flux) divided by the peak flux (Woods & 
Loeb 1994). However, recently the burst lightcurves have become available (Fishman et al. 
1994), and so it has been possible to extract T more directly. Here we define T to be the 
full width at half maximum of the peak, i.e., the time interval from the point when the flux 
rises above half of its maximum value to the point when it falls below half of the maximum 
value. In the case of a burst with many peaks, the value of T was computed only for the 
highest peak. Because 7 m i n depends only weakly on T (7 m i n oc J 1-1 / 4 for spectra extending 
to infinite energies and even weaker for a finite energy cutoff), the uncertainty introduced 
in multiply-peaked bursts is small. The method used for determining and subtracting the 
background is discussed in §3.2. Bursts with data gaps (i.e., where the count rate drops 
suddenly to zero for some length of time due to instrumental error) were excluded from the 
sample. 

Figure 5 shows the results of applying the constraints derived in §2.1, along with equa- 
tions (19)-(21), to the 254 bursts in the second BATSE catalog with Cp ea k > 1 cm -2 sec -1 
and no data gaps. The distributions in the minimum Lorentz factor 7 m i n , the maximum 
baryonic mass Af max , and maximum ambient density n max are shown, in addition to the 
distribution in total 7-ray energy for a spherical source (Ej = 4ttD 2 S, where S is the burst 
fluence). The values of E~ /} Af max , and n max should be multiplied by AO/47T for jets that 
cover a solid angle AO. We have taken the observed test photon energy to be et/(l + z) = 100 
MeV, since this is of the order of the highest photon energy to which BATSE is sensitive. 
We show distributions for the case where the burst spectra extend to infinite energies (Fig. 
5a — a 7 ), and for the case where the spectrum cuts off at 100 MeV (Fig. 6a — a 7 ). Each 
histogram is normalized to unit area, and shows the fractional number of bursts per log- 
arithmic interval in a given quantity x. Thus, each histogram is a discrete plot of xP(x) } 
where P(x) is the probability distribution in x for the burst sample. For spectra extending to 
infinite energies (curve (a) in Fig. 2), one has the leading-order proportionalities Ej oc h~ 2 , 

7min oc h~ 1 / 2 e^ 4: , Mm ax oc /i _3 / 2 e t and n max oc h^l 2 e t Note that since 7 m i n is inde- 
pendent of C pea k to within a factor of order unity, the distribution in 7 m i n simply reflects the 
distribution in peak emission duration T. Also note that only n max is strongly sensitive to 
the test photon energy et; the values of n max would be lowered appreciably for test photon 
energies much higher than 100 MeV. However, a modest increase in 7/7 m i n can compensate 
for this reduction since, to leading order, n max oc 7 s . For the case of an energy cutoff at 100 
MeV (curve (6) in Fig. 2), the dependences on et and h are weaker. The narrow distribution 
for 7 m i n in figure 66 can be understood by inspecting figure 2. One obtains only a small 
range of 7 m i n values by sliding the dashed curve for i? max up and down and changing its 
point of intersection with the lower solid curve for R m m(b). 
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Fig. 5: Probability distributions P(x) for (a) the radiation energy E~ /} (b) 

the minimum Lorentz factor 7 m i n , (c) the maximum baryonic mass Af max , 

and (d) the upper bound on the ambient gas density n max , for the case where 

GRB spectra extend to infinite energies. The histograms show xP(x) for 

the variable x whose logarithm appears on the horizontal axis. We use a 

— 7/4 

test photon energy et = 100(1 + z) MeV. Note that only n max oc e t is 
significantly sensitive to et. 
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Fig. 6: The same as in figure 5, for the case where GRB spectra trun- 
cate above fOO MeV. The dependence of n max on 7 in this case is more 
complicated than in figure 5. 



Table f gives the mean (x), the standard deviation Ax = ((x — (x)) 2 } 1 / 2 , and the log- 
skewness si og = ((log x — (log x)) 3 )/((log x — (log x}) 2 } 3 / 2 for the distributions shown in figures 
5 and 6. A rough estimate of the significance of the skewness is obtained by comparing it to 
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(15/A^tot) 1//2 ~ 0.24, where iVtot = 254 is the total number of data points (Press et al. 1992). 
All of the distributions are negatively skewed in log space, but not very significantly; they 
are well approximated by a log-normal profile. 



Moment 


E 1 l 'erg 


Tmin 


i.W max /.\/. 


^max/cm 3 


Mean 


3.9 x 10 51 


4.6 x 10 2 


0.61 x 10" 6 


2.1 x 10 3 


Std. Dev. 


5.9 x 10 51 


1.4 x 10 2 


1.0 x 10" 6 


2.9 x 10 3 


Log- Skew 


-0.33 


-0.05 


-0.32 


-0.16 



Table la: Mean, standard deviation, and log-skewness of the distributions given in figure 
5, for the case where GRB spectra extend to infinite energies. Values of ^n max are shown 
for 7 = 7 min . 



Moment 


E 1 l 'erg 


Tmin 


i.W max /.\/. 


^max/cm 3 


Mean 


3.9 x 10 51 


8.5 x 10 1 


2.6 x 10" 5 


8.8 x 10 8 


Std. Dev. 


5.9 x 10 51 


0.37 x 10 1 


4.0 x 10" 5 


2.2 x 10 9 


Log- Skew 


-0.33 


-0.05 


-0.33 


-0.20 



Table lb: Mean, standard deviation, and log-skewness of the distributions given in figure 6, 
for the case where GRB spectra cut off at 100 MeV. Values of <^n, max are shown for 7 = 7 m i n . 
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3.2 The Maximum Curvature Radius i? max 

There are now over 500 bursts whose lightcurves are available in the second BATSE 
catalog. The time histories are given in most cases for ~ 240 sec after triggering, although 
some have more (~ 560 sec for some of the later bursts) and some have less. The data 
are summed over the energy channels 50-300 keV, and are also summed over all triggered 
detectors. The time resolution is 64 msec. 

The quantity 8t / \8C /C) in equation (18) may be considered an operational definition 
for the variability timescale of a burst, i.e. the time it takes for the count rate to change by 
an amount comparable to its mean. It is not a unique definition, since it depends on the 
choice of 8t and over which part of the burst lightcurve one places the "time window" of 
width 8t . For our purposes it is in our interest to take 8t as small as possible and choose 
that part of the burst when 8C/C (taken over a time 8t ) is the largest. This gives us the 
strictest possible constraint in equation (18). The smallest value of 8t is limited by the 
time-resolution of the detector. Thus, the constraint on R could be improved substantially 
with better instruments in the future. 

The absolute normalization of the count rates is irrelevant since it is only the ratio 8C / C 

which appears in equation (18). However, to obtain the mean count rate C it is necessary to 

determine and subtract the background count rate. This is difficult to do rigorously because 

of Poisson noise and the fact that there is no systematic way to isolate a stretch of flat 

background for all of the bursts. We therefore estimate the background by binning the data 

very coarsely every ~ 6.4 sec, finding the minimum count rate over the entire burst lightcurve, 

and dividing by 100 to get the DC background count rate Cb ac k f° r 64-msec binning. This 

1/2 

procedure yields an underestimate because of the noise (given by <5Cb ac k ~ ^back)' on ^y 
a very small underestimate since the background is estimated for a very large binning size. 

The noise also poses a problem in determining the maximum value of 8C/C over the 

burst for a given 8t . In the absence of noise, one could simply take the smallest possible 

value of 8t = 64 msec. However, the signal-to-noise ratio is roughly proportional to (<5t ) _1 / 2 , 

so we expect that as we increase 8t , we should improve the purity of the burst signal. We 

nominally require a minimum signal-to-noise ratio of 5. The signal-to-noise ratio is defined 

as S/N = (^C) max /<5Cback? where (<5C) max is the maximum value of 8C for a given 8t and 
1/2 

<*> Cback ~ ^back * s ^ ne n °i se l evei ? given by the square root of the background count rate in 
photons sec -1 . We thus choose the smallest value of 8t such that S/N > 5. If a 8t of 
greater than 0.64 sec is needed, the burst is discarded. The effect of varying the minimum 
S/N is discussed below. 

Figure 7a shows the distribution in values of i? max for the 319 bursts in the second 
BATSE catalog for which S/N > 5 with 8t < 0.64 sec. The value of i? max for each burst 
is calculated with the smallest value of 8t which gives S/N > 5. More than half of the 
bursts have S/N > 5 for 8t = 0.064, at the limit of BATSE's resolution. This, along with 
the fact that 8C/C < 1 by definition, accounts for the "pile-up" of events at low values of 
Rmax in figure 7. Better time resolution will lower the values of i? max , thus strengthening 
our constraints. 

Changing the signal-to-noise cutoff can in principle change the shape of the i? max distri- 
bution. If the cutoff is increased, the sample size will decrease and bursts with larger 8t will 
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Fig. 7: Distribution in values of i? max , the maximum curvature radius 
derived in §4, with a minimum signal-to-noise of (a) S/N = 5 and (6) 
S/N = 10. Values were calculated for 7 = 10 2 , c sg = c, and rj = 1 (cf. Eq. 
(18)). 



constitute a larger part of the sample. This will skew the i? max distribution toward larger 
values. However, the shape of the histogram does not change qualitatively (see Fig. 76), and 
the range of i? max values spanned by the sample does not change. 

The biggest uncertainty in figuring the value of i? max originates from the unknown value 
of the Lorentz factor 7. In principle, one can hope to place an upper bound on 7. Unfortu- 
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nately, no such bound is strict enough to place a useful limit on R. An upper bound on the 
value of 7 achievable in a fireball model is the value reached when the Compton drag time 
exceeds the comoving expansion time in the fireball (Meszaros, Laguna, & Rees 1993). How- 
ever, this constraint depends on the initial radius of the fireball, which is a model-dependent 
quantity. 

In principle, the signal propagation speed be nonrelativistic. If the average GRB 

photon is observed to have an energy ~ MeV, its energy in the rest frame of the expanding 
material is ~ MeV/7. Since 7 ^ 10 2 for most bursts, this would be consistent with a 
nonrelativistic temperature, and hence a nonrelativistic sound speed, in the material rest 
frame. If c sg <C c, then the constraint in equation (18) is strengthened considerably relative 
to equation (11). 



The empirical evidence for a cosmological origin of 7-ray bursts is compelling (Mao 
& Paczyhski 1992; Piran 1992; Meegan et al. 1993; Norris et al. 1993). However, one 
important ingredient is missing: no host galaxies have been found near any well-localized 
bursts (Schaefer 1992; Fenimore et al. 1993). In this section we make the assumption that 
bursts occur in galaxies, and place constraints on the luminosities of the host galaxies and 
on the width of the GRB luminosity function. 

For simplicity, we assume non-evolving populations of both galaxies and 7-ray burst 
sources and a flat (0 = 1) universe. The consequences of introducing galaxy evolution are 
discussed in §5. If 7-ray bursts are associated with neutron stars, there is no reason to expect 
the shape of the GRB luminosity function to evolve with cosmological time, although the 
comoving density of GRB sources may evolve. 

Suppose that the positional error box for a given burst is searched for galaxies in the 
wavelength band A, and that no galaxies are found with an apparent magnitude brighter 
than m. Then, the host galaxy for this burst cannot be brighter than iv max , where 



Here, M\q is the absolute magnitude of a lO 10 /^© galaxy (M10 = —19.52 and —20.69 for the 
B and R bands, respectively), -D(^) is the luminosity distance to the burst source at redshift 
zj, and K\(z) is the /^-correction for a galaxy at redshift z in the wavelength band A. The 
/^-correction simply accounts for the redshifting of photons from shorter wavelengths in the 
galaxy spectrum into the observed band. The correction due to extinction by dust in the 
Milky Way is assumed to be small. The values of z\, are calculated by assuming that bursts 
are standard candles in peak flux and using equations (19) and (20), with the best-fit peak 
emission rate f = 1.0 X 10 57 sec -1 . The peak flux values were obtained from Fenimore et 



Table 2 lists Schaefer's limiting magnitude and band, our standard-candle burst redshift 
zi, the values of Zy max , the size AO of the error box in arcmin 2 on the sky, and the number 
of galaxies iVg a i brighter than the limiting magnitude that one would expect to see in a 
randomly placed field the size of the error box, for eight of the PVO (Pioneer Venus Orbiter) 



4. LIMITS ON THE ASSOCIATION OF GRB WITH GALAXIES 





al. (1993). 
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bursts examined by Schaefer (1992). Values of iv max are given separately for elliptical and 
spiral galaxies, since the /^-correction depends on the galaxy type. The /^-corrections were 
obtained by interpolating between the redshift data points given by Frei and Gunn (1994). 
Our values of zj, are systematically higher than those obtained by Fenimore et al. (1993) by 
Azi, ~ 0.05, since our best-fit luminosity and assumed spectral shape for GRB are different. 
The expected number of galaxies was calculated for the B band from the data of Lilly (1993). 
The numbers for the R band were calculated using a no-evolution model. The sizes of the 
positional error boxes were obtained from the catalog of Atteia et al (1987). Note that the 
nondetection of galaxies near the bursts of 24 Nov 1978 and 6 Apr 1979 is statistically highly 
unlikely, and so we find it surprising that no galaxies were seen. 



Burst 


m (band) 




-^max(ElL) 

(io lo x ) 


-^max(Sp.) 

(io lo x ) 


AO/(l') 2 




19 Nov 78 


16.3 {B) 


0.11 


17.2 


13.3 


8. 


0.01 


24 Nov 78 


20.0 (R) 


0.13 


0.18 


0.17 


48. 


6. 


25 Mar 79 


20.0 (R) 


0.19 


0.42 


0.37 


2 


0.2 


29 Mar 79 


18.0 (B) 


0.25 


31.0 


20.3 


41. 


0.2 


06 Apr 79 


24.7 (B) 


0.17 


0.02 


0.02 


0.26 


3. 


18 Apr 79 


20.5 (B) 


0.14 


0.57 


0.43 


2.9 


0.3 


13 Jun 79 


19.5 (R) 


0.27 


1.51 


1.25 


0.76 


0.05 


16 Nov 79 


21.0 (B) 


0.11 


0.21 


0.17 


3.7 


0.7 



Table 2: Limiting magnitude m (Schaefer 1992), source redshift zj, (calculated using Eq. 
(19)-(20)), limiting luminosity L m?LX in units of 10 W L Q (cf. Eq. (22)), size AO of the 
positional error box in arcmin 2 , and expected number of galaxies iVg a i in the burst error box, 
for eight well-localized PVO bursts. The values of z\, and Zy max are based on the assumption 
that GRB are standard candles in peak flux. Values are shown separately for elliptical (Ell.) 
and spiral (Sp.) galaxies because of the different /^-corrections for each type. 

If GRB do occur in galaxies, one may calculate the probability P of not seeing the 
host galaxy for a given burst in a wavelength band A down to a limiting magnitude m, or 
equivalently, the fraction of the total number of potential host galaxies which are intrinsically 
fainter than the limiting luminosity Zy max given by equation (22). Again we assume no 
evolution and take the luminosity functions for each of the galaxy types to be Schechter 
functions, so that the number galaxies per Mpc 3 in the luminosity interval (L, L + dL) at a 
redshift z is given by 

</>i(z,L) dL = (l + z) 3 (-^-) exp(--^-)^L, (23) 

where i refers to the morphological type and <f>*j, L* jl} and a l for each type taken from 
the fits to the CfA redshift survey data (Marzke, Geller, & Huchra 1994). These fits were 
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obtained using Zwicky magnitudes; we use the same fitted values of <f>*j, etc. for the B } V, 
and R bands probed by Schaefer (1992). The factor of (1 + z) 3 in equation (23) accounts for 
the scaling of the proper density of galaxies with cosmological redshift. Suppose that GRB 
only occur in galaxies brighter than L m [ n . If bursts are standard candles in peak flux, then 
the desired probability is 

yv rL^,i(z b ,m) dL fcU L \ 

where X m ax,i(^6) is given by equation (22) and zj, is the burst redshift. For standard candles, 
zj, can be determined from equations (19) and (20). All of the bursts considered here are 
relatively bright, and have standard-candle redshifts zj, ^ 0.25. This justifies our neglect of 
evolution for the galaxy luminosity function. At such low redshifts, one might also expect a 
low value of P, i.e., most galaxies should be visible down to a limiting magnitude as faint 
as 24. However, bursts are not necessarily standard candles, and the width of the GRB 
luminosity function has not been well-constrained (Lubin & Wijers 1993; Ulmer & Wijers 
1994; Cohen & Piran 1995). In principle, the probability that apparently bright bursts 
originate at large distances increases as the width of the GRB luminosity function increases. 
It is therefore conceivable that above a certain width of the GRB luminosity function, the 
probability for non-detection of the distant host galaxies would be considerable. 

To quantify the width of the GRB luminosity function over a large range of values we 
use a "log-normal" distribution. The probability that a given burst falls in the interval 
(T, T + dT) is then 

where V is the intrinsic peak emission rate in photons sec 1 and To is the best-fit value of V 
obtained in §3.1 (cf. Eq. (19)). The value of Pb(T) in equation (25) is normalized to a unit 
area. A burst with observed peak flux Cp ea k then has a probability pb(z)dz = pb(T)(dT / dz)dz 
of being in the redshift range (z,z + dz). For an = 1 cosmology, the redshift probability 
distribution is thus 

( \ ^Peak e-" 2/2 n , 3 r— 1 f -[ln(T(z) /Tp)] 2 \ 

Pb( z ) = —F, ^ (1 + z ~ o V 1 + z + -) exp <^ — \ , (26) 



Co 2 V 2' r [ 2a 2 



where Y(z) = 167r(c/# ) 2 C peak (l + z - vTT*) 2 (cf. Eq. (19) and (20)), and C = 
r /[167r(c/# ) 2 ]. In this case, the probability in equation (24) is modified to an integral 
over the burst redshift distribution, 



z 



Ei JIT dz Ph( z ) dL <M Z ' z ) 

Figure 8 shows the probability P that the host galaxies were not seen for six of the 
eight PVO bursts under consideration. The burst luminosity function (25) was chosen to be 
centered on the best-fit emission rate To from §3, with a width a ranging between 0.1 and 



22 




Fig. 8: Probability P of not seeing the host galaxy as a function of the 
width a of the GRB luminosity function (cf. Eq. (25)), for six of the PVO 
bursts considered by Schaefer (1992). Values are plotted for L m [ n = 10 8 Lq 
(dotted line), lO 9 /© (dashed line), and 10 10 Lq (solid line). The bursts of 
19 Nov 1978 and 29 Mar 1979 both have P 1 for all values of a (see Table 
2)- 



2.0. Note that this range is larger than that given by good fits to the N vs. Cp ea k curve 
in Fig. 4c. The morphological classes for which we obtained our /^-corrections from Frei 
& Gunn (1994) were assumed to correspond roughly to the classifications used by Marzke 
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et al. (1994) in determining the galaxy luminosity functions. The probabilities are shown 
as a function of a for each burst, for the three values L m [ n = lO 8 iv0, lO 9 iv0, and lO lo iv0. 
The events of 19 Nov 1978 and 29 Mar 1979 both have P 1 even for a = (i.e., standard 
candles), and are therefore not shown. For most of the other bursts, the probability becomes 
large only for values of a which are inconsistent with the fits in figure 4c. We expect that 
choosing a GRB luminosity function which gives a better N vs. C pea k fit for large widths 
will not be sufficiently skewed toward high luminosities to produce a stronger dependence of 
P on the luminosity function width. The fact that no galaxies were seen in the error boxes 
for the bursts of 24 Nov 1978 and 6 Apr 1979 is surprising even if one does not assume an 
association between GRB sources and galaxies. This leads us to question the validity of 
some of the observed limits. 

5. SUMMARY AND DISCUSSION 

In this paper, we have derived several empirical constraints on the physical properties 
of 7-ray burst (GRB) sources. In §2.1 and §3.1, the requirement of optical thinness to pair 
production was used to constrain the bulk Lorentz factors and other physical properties of 
bursts in the second BATSE catalog. In §2.2 and §3.2, the time variabilities of bursts in the 
second BATSE catalog were used to place an upper bound on the curvature radius of the 
emitting surface which produces the burst. In §4, we used the data on several PVO bursts 
(Schaefer 1992; Fenimore et al. 1993) to examine the likelihood that GRB occur in unseen 
galaxies as a function of the width of the GRB luminosity distribution. 

The results in §2.1 and §3.1 were obtained under four simplifying assumptions. First, 
GRB were taken to be standard candles in peak luminosity, with C pea k as a distance indicator. 
In fact, the K-S significance of the N — C pea k fit (Fig. 4c) is not very sensitive to the width of 
the GRB luminosity function (Lubin & Wijers 1993; Ulmer & Wijers 1994; Cohen & Piran 
1995). As discussed in §4 and below, a wide GRB luminosity function may be needed to 
explain the nondetection of galactic hosts near burst locations on the sky. In this case, a 
single value of C pea k would correspond to a distribution in redshifts, and each burst would be 
spread over more than one bin in the histograms of figures 5 and 6. The second assumption 
was that the peak flux emission of the source is steady over a timescale ~ Rq / c. A shorter 
burst would have a smaller optical depth, because photons emitted at an angle ~ I/7 would 
be unable to catch up with photons emitted along the line of sight to the center of the source. 
If this assumption is violated, the constraint on 7 is weakened. Finally, the geometric time 
delay expressed in equation (11) was derived assuming that the emission comes from a 
surface of size -R0/7. If the emitting region is a smaller clump of matter, the constraint on 
7 is weakened. 

Bearing these assumptions in mind, we arrive at a few physically interesting results from 
§2.1 and §3.1. The energy distribution in figures 5a and 6a ends at about the binding energy 
of a neutron star, ~ 10 53 ergs. This result is consistent with a variety of models that associate 
cosmological 7-ray bursts with neutron stars (Eichler et al. 1989; Meszaros & Rees 1992; 
Narayan, Paczyhski, & Piran 1992; Usov 1992; Loeb 1993), although in some models it is 
difficult to extract all the binding energy in 7-rays. The distribution of maximum ambient 
gas densities n max (Figs. 5a 7 and 6a 7 ) allows bursts to reside within the interstellar medium 
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of galaxies, where the density n ISM 1 cm , and is also consistent with a molecular cloud 
environment (n ~ 10 2-4 cm -3 ). 

The limit i? max derived in §2.2 and §3.2 on the curvature radius of the emitting surface 
was obtained assuming that the causally disconnected regions on the surface act as truly 
independent sites. It is possible that many such regions are lit up at the same time due 
to some overlap in their past light cones, but we ignore this possibility, assuming that the 
conditions responsible for the emission are local. The fireball model of Meszaros, Laguna, & 
Rees (1993) places an effective lower bound on the value of R at which the main emission 
can take place. Since the observed 7-ray energy must be converted from the bulk kinetic 
energy of the fireball matter, the main part of the burst should not occur until the expansion 
of the fireball is decelerated appreciably by the inertia of the surrounding material. This 
happens when the mass of interstellar matter swept up by the fireball is comparable to the 
bulk kinetic energy of the baryons in the fireball. For a spherical fireball with a typical 
energy (~ 10 51 ergs) and interstellar ambient density (~ 1 cm -3 ), the radius at which the 
deceleration becomes significant is ~ 10 16 ( 7 /10 3 ) -2 / 3 cm. To reconcile this with the typical 
value of i? max = 10 14 ( 7 /10 2 ) 2 found in figure 7, one needs 7 ^ 10 3 , a stricter bound than 
that derived from the optical depth to pair production in §3.1. The alternative is to have 
emission from collisions between shells in the fireball at smaller radii (Rees & Meszaros 
1994). If the signal speed in the local rest frame is much smaller than the speed of light, the 
constraints on i? max and 7 become stronger. It should be emphasized that i? max applies to 
Galactic as well as cosmological sources. 

In §4 it was assumed that the luminosity functions of galaxies did not vary with redshift 
(the "no-evolution" model). The number of galaxies actually observed at blue magnitudes 
> 20 is larger than the no-evolution prediction. Thus, introducing an evolutionary model 
which agrees with galaxy counts will only strengthen our result that the probability of not 
seeing the GRB hosts is low. The corrections due to extinction by dust in the Milky Way 
were assumed to be small. This is justified for five of the eight bursts considered since these 
bursts occurred at galactic latitudes \b\ > 55°, where there is virtually no extinction; for the 
other three bursts, the extinction correction is not likely to be much more than two-tenths of 
a magnitude (Burstein & Heiles 1982). The non-detection of galaxies near the bursts of 24 
Nov 1978 and 6 Apr 1979 is surprising, given the expected counts based on galaxy surveys. 
We would expect to find a few galaxies brighter than the limiting magnitude for each of 
these two burts in a randomly-placed field the size of the associated error box (Table 2). In 
the most restrictive case, the event of 6 Apr 1979, the host galaxy cannot be brighter than 
~ lO 8 iv0 if GRBs are standard candles (Table 2). Even if one allows for a GRB luminosity 
function centered on the best-fit value from §3.1 with a width as large as its mean, the 
probability of not seeing the host galaxy for this burst is very small unless galaxies with 
L ^ lO 8 iv0 are more likely to host 7-ray bursts. Generally, one might expect that the 7- 
ray burst rate in a given galaxy would be roughly proportional to the number of stars and 
hence to the total luminosity. However, if most of the light in the universe originates from 
faint galaxies, one should see many GRB events occurring in faint host galaxies. In our 
no-evolution model, the faint-end logarithmic slope of the total galaxy luminosity function is 
a ~ — 1 (Marzke et al. 1994), so that the luminosity per logarithmic interval is L 2 (f>(L) oc L, 
and thus faint galaxies do not contribute most of the light. However, an excess of faint 
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galaxies at redshifts z ^ 0.1 or a variation in the mass function of stars in faint galaxies 
may allow for a larger abundance of neutron stars in faint hosts. In addition, there may 
be a large population of sources which exist outside of galaxies. High- velocity (<; 500 km 
sec -1 ) neutron star binaries which escape from galaxies could easily wander outside of the 
1 arcmin 2 box within a Hubble time. There could also be a population of stars formed in 
intergalactic space. As shown in figure 8, another solution to the problem is to introduce a 
broad GRB luminosity function. However, figure 4c demonstrates that any such luminosity 
function is significantly constrained by the number-count statistics of the bursts. 

We thank Stirling Colgate, George Field, Josh Grindlay, Shude Mao, Ramesh Narayan, 
Bohdan Paczyhski, and Eli Waxman for insightful comments and suggestions. 
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